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Abstract. We present the redshift evolution of the restframe galaxy luminosity function (LP) in th e red r', i', and 
z' band s as derived from the FORS Deep Field (PDF), thus extending the results published in iGabasch et all 
to longer wavelengths. Using the deep and homogeneous I-band selected dataset of the PDF we are able 
to follow the red LPs over the redshift range 0.5 < z < 3.5. The results are based on photometric redshifts for 
5558 galaxies derived from the photometry in 9 filters achieving an accuracy of Az/{zspec + 1) < 0.03 with only 
~ 1% outliers. A comparison with results from the literature shows the reliability of the derived LPs. Because of 
the depth of the PDF we can give relatively tight constraints on the faint-end slope a of the LP: The faint-end of 
the red LPs does not show a large redshift evolution and is compatible within la to 2a with a constant slope over 
the redshift range 0.5 < ^ < 2.0. Moreover, the slopes in r', i', and z' are very similar with a best fitting value 
of a = —1.33 ± 0.03 for the combined bands. There is a clear trend of a to steepen with increasing wavelength: 
ctuvt^u' = —1.07 ±0.04 — > Og/fes = —1.25 ±0.03 — > Qr'&i'&z' = —1.33 ±0.03. We subdivide our galaxy sample 
into four SED types and determine the contribution of a typical SED type to the overall LP. We show that the 
wavelength dependence of the LP slope can be explained by the relative contribution of different SED-type LPs 
to the overall LP, as different SED types dominate the LP in the blue and red bands. Furthermore we also derive 
and analyze the luminosity density evolution of the different SED types up to z ~ 2. 

We investigate the evolution of M* and </>* by means of the redshift parametrization M* (z) — Mg + a ln(l ± z) and 
<j)*{z) = (I>^{1 + ^)^ Based on the PDF data, we find only a mild brightening of M* (a^' ~ —0.8, and a^/ ^/ ~ —0.4) 
and decrease of 0* (6^/ ^/ ^/ ~ —0.6) with increasing redshift. Therefore, from (z) ~ 0.5 to (z) ~ 3 the characteristic 
luminosity increases by ~ 0.8, ~ 0.4 and ~ 0.4 magnitudes in the r', i', and z' bands, respectively. Simultaneously 
the characteristic density decreases by abo ut 40% in all ana lyzed wavebands. A comparison of the LPs with 
semi-analytical galaxy formation models bv lKauffmann et alJ lfl99ftl shows a similar result as in the blue bands: 
the semi-analytical models predict LPs which describe the data at low redshift very well, but show growing 
disagreement with increasing redshifts. 
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1. Introduction 

One of the major efforts in extragalactic astronomy is to 
derive and analyze the restframe galaxy luminosity func- 
tion in different bandpasses and redshift slices in order to 
follow the time evolution of the galaxy populations by a 
statistical approach. This is of particular importance be- 
cause the energy output at different wavelengths is domi- 
nated by stars of different masses. While galaxy luminosi- 
ties measured in the ultraviolet are sensitive to the energy 
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output of hot, short-living O and B type stars and there- 
fore to ongo ing star formation l)Tinslevll97lHMadau et all 
Il996lll998|) . the optical and NIR luminositie s provide con- 
strain ts on more evolved stellar populations ijHunter et alJ 

This can be used, in principle, to derive the evolu- 
tion of b asic galaxy propert ies as the stellar mass function 
(see e.g. iDrorv et al.ll2005l and re ferences therein), the 
star f ormation rate density (see e.g. lPerez-Gonzalez et all 

and re ferences therein) or t he specific star formation 
rate (see e.g. lFeulner et al.ll2005[ and references therein). 
The determination of these quantities, however, is based 
on assumptions, e.g. on the shape of the initial mass func- 
tion or on the details in modeling the stellar population 
like age, chemical composition, and star formation history. 
Hence studying the LF at different wavelengths and cos- 
mic epochs offers a more direct approach to the problem 
of galaxy evolution. 

As the LF is one of the fundamental observational tools, 
the amount of work spent by different groups in deriv- 
ing accurate LFs is substantial. Based on either spec- 
troscopic redshifts, drop-out techniques, or photometric 
redshifts, it has been possible to derive luminosity func- 
tions at differ e nt red s hifts in the ultravio le t fc blue band 
('Baldrv et aP 120 05': 'Croton et aP l2005t lArnouts et alJ 
2005: Biidavari et al. .200,^ .Trever et a] . I [200,51 see also 



Gabasch'et^n2004al and r e ferenc e s therein), i n the 



red bands | Lin et al.' 'l996', 'l997; Brow n et al 
Shaplev et al.> i2001^ .Wolf et al.> .2003.: Chen et al 



Ilbert et all l2004t iDahlen et all 120051 iTrentham et all 
20051) as well as in the near -IR bands ("LQ^g^^^l^^^^ ; 



2001 



20031 



Kochanek et a,].l200lHGo]e et aL.2001: Balogh et alJ2001 



Drorv et abl^ lHua,ng et a1.ll2003HFeu]ner et a1.ll2003l 

i 



, , Q£ 

Pozzetti et 311120031 iDahlen et aDl2005|) . 

The evolution of the characteristic luminosity and den- 
sity of galaxy popula tions can be an alyzed by fitting 
a Schechter hmction l|Schechterl Il97fil) to the LF. The 
redshift evolution of the three free parameters of the 
Schechter function, the characteristic magnitude M*, the 
density 0*, and the faint-end slope a can be used to 
quantitatively describe the change of the LF as a func- 
tion of redshift. Unfortunately, the Schechter parametriza- 
tion of the LF cannot account for possible excesses at 
the bright and faint end or other subtle shape deviations. 
Furthermore, the Schechter parameters are highly corre- 
lated making it challenging, but not impossible, to clearly 
separate the ev olution of the different parameters (see e.g. 
lAndreonlliooi for a discussion). 

The evolution of the LFs is also very suitable to con- 
strain the free parameters of theoretical models (e.g. semi- 
analytical or smoothed particle hydrodynamics models). 
Ideally a comparison between model predictions and ob- 
servations should be done simultaneously for different 
wavebands (UV, optical, NIR) and for different redshift 
slices as different stellar populations are involved in gener- 
ating the flux in t he different bands. Therefore, the FDF 
llHeidt et alJ 120031) provides a unique testing ground for 
model predictions, as the depth and the covered area al- 
low relatively precise LF measurements from the UV to 



the z'-band up to high redshift in a very homogeneous 
way. 

In this paper we extend the me asurements of the blue 
luminosity functions presented in iGabasch etall l|2004al 
hereafter FDFLF I) to the red r', i', and z' bands. In Sect.H 
we derive the LFs and show the best fitting Schechter pa- 
rameters M*, (/)*, and a in the redshift range 0.5 < z < 3.5. 
We also present a detailed analysis of the slope of the LF 
as a function of redshift and wavelength. Furthermore, we 
analyze the contributions of different SED types to the 
overall LF and present the evolution of the type depen- 
dent luminosity density up to redshift 2 ~ 2. Sect. |31 shows 
a parametric analysis of the redshift evolution of the LF, 
whereas a comparison with the LFs of other surveys as 
well as with model predictions is given in Sect. 01and in 
Sect. 121 respectively. We summarize our work in Sect. El 
Throughout this paper we use AB magnitudes and 
adopt a A cosmology with fl^j = 0.3, JIa = 0.7, and 
Hq = 70kms-iMpc~^ 

2. Luminosity functions in the r', i', and z' bands 

The results presented in this paper are all based on the 
deep part of the I-band selected catalog of the FDF 
l(Heidt et alJl2003l) as introduced in FDFLF I. Galaxy dis- 
tances are determined by t he photometric redshift tech- 
nique ijBender et al.l l200l(l with a typical accuracy of 
Az/(zs„ pr. -I- 1) ^ 0.03 if compared to the spectroscopic 
sample (jNoll et al.ll2004l iBohm et al1l2004() of more than 
350 objects. To derive the absolute magnitude for a given 
band (which will be briefly summarized below) we use the 
best fitting SED as determined by the photometric red- 
shift code, thus reducing possible systematic errors which 
could be introduced by using k-corrections applied to a 
single observed magnitude. To account for the fact that 
some fainter galaxies are not visible in the w hole survey 
volume we perform a V/Vmax l|Schmidtlll968l) correction. 
The errors of the LFs are calculated by means of Monte- 
Carlo simulations and include the photometric redshift er- 
ror of every single galaxy as well as the statistical error 
(Poissonian error). To derive precise Schechter parame- 
ters we limit our analysis of the LF to the magnitude bin 
where V/Vmax < 3. We also show the uncorrected LF in 
the various plots as open circles. We do not assume any 
evolution of the galaxies within the single redshift bins, 
since the number of galaxies and the distance determina- 
tion based on photometric redshifts would not be able to 
constrain it. The redshift binning was chosen such that we 
have good statistics in the various redshift bins and that 
the influence of redshift clustering was minimized. In or- 
der to have good statistics at the bright end (rare objects) 
of the LF we had to slightly change some of the redshift 
bins if compared to FDFLF I. The new redshift binning 
together with the number of galaxies in every bin is shown 
in Table As can be seen from Table ^ the redshift in- 
tervals are approximately the same size in ln(l -I- z) and 
most of the results we are going to discuss in this paper 
are based on 700 - 1000 galaxies per redshift bin. 
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Table 2. Slope of the LF for all wavelengths and all redshifts as derived from a 3-parameter Schechter fit. 



z 


a (r') 


a(i') 


a (z') 


0.45 - 0.85 
0.85 - 1.31 
1.31 - 1.91 
1.91 - 2.61 
2.61 - 3.81 


-1.37 (+0.04 -0.04) 
-1.25 (+0.06 -0.04) 
-1.30 (+0.16 -0.09) 
-1.01 (+0.15 -0.14) 
-0.98 (+0.17 -0.17) 


-1.37 (+0.04 -0.03) 
-1.27 (+0.06 -0.05) 
-1.50 (+0.13 -0.10) 
-1.03 (+0.17 -0.14) 
-1.03 (+0.15 -0.13) 


-1.39 (+0.04 -0.04) 
-1.34 (+0.06 -0.04) 
-1.45 (+0.12 -0.09) 
-0.97 (+0.17 -0.12) 
-1.01 (+0.15 -0.13) 



Table 1. Number of galaxies in the FDF for the redshift 
intervals used for computing the LFs. Note that we de- 
rive the LF in all redshift bins, but exclude the lowest 
(z < 0.45) and highest redshift bin (z > 3.81) from our 
analysis of the LF evolution, since the lowest redshift bin 
corresponds to too small a volume while the z > 3.81 bin 
suffers from extrapolation errors. 



redshift 


number 


fraction 


interval 


of galaxies 


of galaxies 


0.00 - 0.45 


808 


14.54 % 


0.45 - 0.85 


1109 


19.95 % 


0.85 - 1.31 


1029 


18.51 % 


1.31 - 1.91 


880 


15.83 % 


1.91 - 2.61 


816 


14.68 % 


2.61 - 3.81 


718 


12.92 % 


> 3.81 


196 


3.53 % 


unknown 


2 


0.04 % 



2.1. The slope of the LF as a function of redshift 



depth). Under the assumption that a does not depend 
on redshift. Table O (upper part) yields the slopes' 
best error-weighted values in the redshift range from 
{z) ^ 0.65 to (z) ^ 1.6 (including also the higher redshift 
bins changes a only marginally). Since the slopes in all 
bands are very similar we derive a combined slope of 
Or'&i'&z' = —1.33 + 0.03 (Table|21 lower part). 
Almost all of the slopes listed in Tableware compatible 
within la — 2a with a = — 1.33±0.03. Therefore, we fixed 
the slope to this value for the further analysis. Please 
note, that this slope is steeper than for the blue bands 
{auvScu' = —1.07 and Ug'^B = —1-25), but it follows the 
trend observed in FDFLF I: With increasing wavelength 
the slope steepens, i.e. the ratio of faint to bright galaxies 
increases. This trend is illustrated best in Fig. ^ where 
we combine the results derived in FDFLF I with those 
of this work and plot the wavelength dependence of 
the LF slope. As we will show in Sect. 12.31 this effect 
can be explained by the contribution of different galaxy 
populations to the overall LF in the various wavebands. 



Table 3. In the upper part of the Table we show the slope 
a of the luminosity functions for the different wavebands 
as determined from an error- weighted fit to the data under 
the assumption that a{z) = const. In the lower part of the 
Table we show the best values of a after combining the 
data of all bands. 



filter 


ct{z) — const. 


r' 


-1.30 + 0.05 


i' 


-1.33 + 0.05 


z' 


-1.35 + 0.05 


r' & i' & z' 


-1.33 + 0.03 



2.2. I selection versus l+B selection 

We checked the dependence of our results on the selec- 
tion band by comparing the I band selected catalog and 
the I+B selected FDF catalog . The combined catalog has 
been described in iHeidt et alJ l|200.ltl and reaches limiting 
magnitudes of / ~ 26.8 and B ~ 27.6. In the combined 
sample M* agrees within its la errors with the values de- 
rived from the I-band catalog only. The slope a tends to 
be slightly steeper in the combined sample, but by not 
more than la. The larger number of objects in the com- 
bined catalog mostly influences the characteristic density 
4>* which is a factor of 1.05 to 1.20 larger (depending on 
the redshift bin). Given the errors of (f>*, this is in the 
order of la to 2a. 



To investigate the redshift evolution of the faint-end 
slope of the LF, we fit a three parameter Schechter 
function (M*, (p* , and a) to the data points. The best 
fitting slope a and the corresponding la errors for the 3 
wavebands are reported in Table |21 for the various redshift 
bins. 

It can be inferred from Table|2]that there is only marginal 
evidence for a change of a with redshift (at least up to 
z ~ 2 where we are able to sample the LF to a suitable 



2.3. The slope of the LF as a function of wavelength 

To better understand the filter-dependence of the LF slope 
shown in Fig. ^ we analyze the contribution of different 
galaxy types to the overall LF. Thus, we subdivide our 
galaxy sample into four SED types and analyze the type- 
dependent LF, i.e. we determine the contribution of a typ- 
ical SED type to the overall LF. The SEDs are mainly 
grouped according to the UV-K color (sec Fig. EJ: for 
increasing spectral type (SED type 1 — > SED type 4) the 
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Fig. 3. LFs for the four SED types in the redshift intervals 0.45 < z < 0.85 (upper panels), 0.85 < z < 1.31 (middle 
panels), and 1.31 < z < 1.91 (lower panels): SED type increases from left to right. The filters are color coded and 
denoted in the upper part of the various panels. For clarity, the three parameter Schechter function fit to the data is 
shown as solid line with the same color coding as the LF. 



SEDs become bluer, i.e. the UV flux (and thus the recent 
star formation rate) increases if compared to the K-band 
flux. IPannella et al.l (jiflQ^ analyzed the morphology of 
about 1400 galaxies in the FDF down to / '--^ 25 mag on 
HST (ACS) data, and find a good correlation between the 
four main SED types and the morphology of the galaxies 
(at least up to redshift z ^ 1.5). 

The four SED types also show a sequence in the restframe 
U-V color often use d to discriminate be tween blue and 
red galaxies (see e.g. iGialloneo et al .120051 and references 
therein). As the restframe U-V color includes the 4000 A 
break it is quite sensitive to galaxy properties as age and 
star formation. The U-V color lies in the range between 
2.3 - 1.9, 2.0 - 1.6, 1.6 - 0.9, and 0.9 - for SED type 1, 2, 
3, and 4, respectively. Therefore, in a rough classification 
one can refer to SED types 1 and 2 (SED type 3 and 4) as 
red (blue) galaxies. We use the same SED cuts at all red- 



shifts (see below), i.e. we do not us e the time evolution o f 
the galaxy color bimodality (see e.g. IGialloneo et al .l2005^ 
to redefine the main SED type of a galaxy as a function 
of redshift. 

We show in Fig.|21the LFs for the four SED types in 
three redshift intervals: 0.45 <z< 0.85, 0.85 < z < 1.31, 
and 1.31 < z < 1.91. The SED type increases from the left 
panel to the right panel, i.e. the extremely star-forming 
galaxies are shown in the rightmost panel. The LFs for 
the different filters are color coded and denoted in the 
upper part of the various panels. We show every LF to 
the limiting magnitude where the V/Vmax begins to con- 
tribute by at most a factor of 1.5, being more conservative 
as for the overall LF {V/Vmax < 3 for every bin). For clar- 
ity, a Schechter function fit to the data is shown as solid 
line using the same color coding as for the LF. 
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Fig. 4. Schechter functions fitted to the LFs in the UV (left panels) as well as in the i'-band (right panels) for the 
redshift intervals 0.45 < z < 0.85 (upper panels), 0.85 < z < 1.31 (middle panels), and 1.31 < z < 1.91 (lower panels). 
The solid lines show the best fitting Schechter functions for the four SED types. The SED type is color coded and 
denoted in the upper left panel. The dotted black line shows the total LF (as fitted to the data) whereas the shaded 
region represent the corresponding la error of the latter. For the slope values of the different SED types see Table El 
Please note that we show the SED type LFs to the limiting magnitude where the V/Vmax begins to contribute by at 
most a factor of 1.5, being more conservative than for the overall LF for which we allow a correction factor of 3. 



First of all it is clear from Fig. 13 that the faint-end of 
the LF is always dominated by SED type 4 galaxies. This 
is true for all analyzed bands. If we focus on the bright 
end of the SED type 4 LFs, we only see a relatively small 
variation between the different filters. On the other hand, 
the difference between the filters for SED type 1 (for the 
bright end) is very large. Although (because of the low 
number density) SED type 1 does not contribute at all 
to the faint-end of the LFs, the picture changes for the 
bright end. While for the bright end of the LF in the UV 



(black line), SED type 1 and 4 galaxies have about the 
same number density, in the red bands SED type 1 galax- 
ies dominate the LF. 

This trend applies for all three redshift bins, although it is 
more pronounced at lower redshift. It explains naturally 
the change of the LF slope as a function of waveband. 
This can be best seen in Fig. 01 where we concentrate 
on only two filters. There we show the Schechter func- 
tions fitted to the LFs in the UV as well as in the i'-band 
for the redshift intervals 0.45 < z < 0.85, 0.85 < z < 1.31, 
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Fig. 5. Luminosity densities in the UV (left panel) and i' (right panel) bands for the four SED types in the redshift 
intervals 0.45 < z < 0.85, 0.85 < z < 1.31, and 1.31 < z < 1.91. The luminosity densities are completeness corrected 
to zero luminosity (ZGL, see text for details). The values are listed in Table01 



Table 4. Luminosity densities in the UV and i' bands for the four SED types. See also Fig. |S1 



niter 


orjU type 


redshift 


luminosity density 


error 


completeness correction (ZGL) 








W Hz"^ Mpc"^ 




% 


UV (2800A) 


1 


0.45 - 0.85 


l.Ole+18 


1.5e+17 


0.1 






0.85 - 1.31 


6.16e+17 


l.Oe+17 


0.7 






1.31 - 1.91 


2.68e+17 


7.6e+16 


8.9 




2 


0.45 - 0.85 


7.58e+17 


1.3e+17 


0.5 






0.85 - 1.31 


6.00e+17 


1.2e+17 


2.0 






1.31 - 1.91 


2.26e+17 


7.0e+16 


9.0 




3 


0.45 - 0.85 


6.88e+18 


7.9e+17 


0.6 






0.85 - 1.31 


6.86e+18 


5.4e+17 


2.0 






1.31 - 1.91 


4.27e+18 


5.6e+17 


4.7 




4 


0.45 - 0.85 


8.29e+18 


5.1e+17 


8.4 






0.85 - 1.31 


1.30e+19 


8.4e+17 


6.9 






1.31 - 1.91 


1.83e+19 


1.9e+18 


21.4 


i' 


1 


0.45 - 0.85 


5.72e+19 


8.7e+18 


0.3 






0.85 - 1.31 


3.22e+19 


5.4e+18 


1.5 






1.31 - 1.91 


1.29e+19 


3.5e+18 


10.2 




2 


0.45 - 0.85 


2.24e+19 


3.9e+18 


0.8 






0.85 - 1.31 


1.61e+19 


3.2e+18 


3.6 






1.31 - 1.91 


6.60e+18 


1.6e+18 


11.6 




3 


0.45 - 0.85 


6.21e+19 


6.9e+18 


0.4 






0.85 - 1.31 


6.56e+19 


5.1e+18 


1.6 






1.31 - 1.91 


3.79e+19 


5.4e+18 


7.4 




4 


0.45 - 0.85 


2.59e+19 


1.5e+18 


5.8 






0.85 - 1.31 


3.84e+19 


2.4e+18 


5.7 






1.31 - 1.91 


5.17e+19 


4.5e+18 


18.9 
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Table 5. Slope of the UV (2800 A) and i' band LF for the different SED types from a 3-parameter Schechter fit. 
The Schechter functions are shown in Fig. 0] 



z 


filter 


Q for SED type 1 


a for SED type 2 


a for SED type 3 


a for SED type 4 


0.45 - 0.85 


UV 


-1.06 (+0.16 -0.10) 


-1.27 (+0.08 -0.04) 


-1.12 (+0.11 -0.07) 


-1.19 (+0.13 -0.11) 


0.45 - 0.85 


i' 


-1.11 (+0.15 -0.02) 


-1.17 (+0.22 -0.03) 


-1.12 (+0.12 -0.11) 


-1.23 (+0.07 -0.10) 


0.85 - 1.31 


UV 


+0.38 (+0.60 -0.37) 


-0.71 (+0.62 -0.27) 


-0.68 (+0.17 -0.15) 


-1.14 (+0.12 -0.08) 


0.85 - 1.31 


i' 


+1.04 (+0.65 -0.68) 


-0.62 (+0.79 -0.32) 


-0.84 (+0.15 -0.13) 


-1.09 (+0.11 -0.06) 



1.4 



1.2 



-1 - 




2000 4000 6000 8000 
wave length [A] 




10* 1C00 

X [A] 



Fig. 1. Slope of the LF as a function of wavelength. The 
filled red squares denote the values derived in FDFLF I 
whereas the filled red dots are taken from this work 
(Table 13 upper part). Local slope valu es (black) are 
shown as filled squ ares llBaldrv et al.ll2005^■ open squa res 
l)Driver et alJ2005l). filled triaiigles llBudavari et alJ20 05). 
open triangles llBlanton et al 
i Norberg et al.ll2002l) . 



il 



and as an asterisk 



and 1.31 < z < 1.91. We plot the single Schechter func- 
tions for all four SED types as well as for the overall LF. 
In the UV the overall LF (dotted line) is completely dom- 
inated by the SED type 4 galaxies. On the other hand 
the overall LF in the i'-band is mainly dominated by SED 
type 1 to type 3 at the bright end, and SED type 4 at the 
faint-end. This results in a steeper slope for the overall 
LF. 

Please note that in Fig. Eland Fig.0]we show the SED type 
LFs and Schechter functions to the limiting magnitude 
where the V/Vmax begins to contribute by at most a fac- 
tor of 1.5, being more conservative than for the overall LF 
for which we allow a correction factor of 3. Furthermore, 
all Schechter functions in Fig.^are fits to the data points. 



Fig. 2. SEDs grouped according to their spectral type. 
See text for details. 

This is also true for the overall Schechter function which 
is not the sum of the individual SED type Schechter func- 
tions and explains why, at the bright end, the overall 
Schechter function is in some plots below individual SED 
type Schechter functions. 

Another interesting aspect which can be inferred from 
Fig. 21 is the fast decrease in number density of bright 
SED type 1 to 3 galaxies if compared to SED type 4 
galaxies (for increasing redshift). Therefore at high red- 
shift (z ^ 2) SED type 4 galaxies start to dominate also 
the overall i'-band LF. This can be seen best, if one fol- 
lows the redshift evolution of the type dependent luminos- 
ity density (LD), i.e. the integrated light emitted by the 
different SED types. The results (for the UV and i' bands) 
are shown in F i g. El ^y e calculated the LD as described in 
ICabasch et alJ 1 2n04bl) : First, we derive the LD at a given 
redshift by summing the completeness corrected (using a 
y/Vmax correction) luminosity of every single galaxy up to 
the absolute magnitude limits. Second, we apply a further 
correction (to zero galaxy luminosity) ZGL, to take into 
account the missing contribution to the LD of the fainter 
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galaxies. To this end we use the best-fitting Schechter 
function for a slope a constant with redshift. For every 
SED type we derive a{z) = const, by an error-weighted 
averaging of the slopes given in Table |S1 This results in 
slopes between a = —0.98 and a = —1.25. For the FDF 
the ZGL corrections are at most 22% in size (see last col- 
umn in TableQ. The small ZGL correction employed here 
stems from the faint magnitude limits probed by our deep 
FDF data set and the relatively flat slopes of the Schechter 
function. Errors are computed from Monte Carlo simula- 
tions that take into account the probability distributions 
of the photometric redshifts and the Poissonian error. 
As shown in the left panel of Fig. [S] the contribution of 
type 1 and 2 galaxies to the UV LD is negligible at all an- 
alyzed redshifts. SED type 3 and 4 completely dominate 
the UV output and although the number density of these 
galaxies decreases with increasing redshift the luminosity 
density (and thus the SFR) increases. 
If we analyze the i'-band LD, in the lowest redshift bin 
SED type 1 and 3 dominate (by a factor of about three if 
compared to type 2 and 4) and have about the same LD. 
At higher redshifts the relative contribution of the differ- 
ent SED types changes because the LD of type 1 and 2 
galaxies decreases with increasing redshift and SED type 
3 and 4 take over. 

A detailed analysis of the type dependent LF will be 
presented in a future paper (Gabasch et. al., in prepa- 
ration) where we combine the I-band selected MUNICS 
catalog (MUNICS IX, Feulner et. al., in preparation, ^ 
900 arcmin^) with the FDF {'^ 40 arcmin^) catalog. This 
overcomes the small volume of the FDF at lower redshift 
making it possible to include also rare bright objects in 
the analysis of the LF. First results in the MUNICS fields 
will be presented in MUNICS IX. 



2.4. The redshift evolution of the LFs 

In this section we analyze the LF by means of a Schechter 
function fit with a fixed slope of a = —1.33. In Fig. 
and Fig. [3 we present the LFs in the r'-band and in the 
i'-band, while the results for the z'-band can be found 
in Fig. |S1 The filled (open) symbols denote the LF with 
(without) completeness correction. The solid lines show 
the Schechter function fitted to the luminosity function. 
The best fitting Schechter parameter, the redshift binning 
as well as the reduced are also listed in each figure. 
The values of the reduced ^-re very good for all red- 
shift bins below z ^ 2. We do not fit our lowest redshift 
bin data {{z) ^ 0.3) with a Schechter function, because the 
volume is too small. For co mparison we also show the local 
LF derived bv lBlanton eTal. (2003) in the SDSS (see also 
Fig.EJ. The best fitting Schechter parameters and corre- 
sponding la errors are summarized in Table El Table [T] 
and Table|Slfor the r', i', and z' bands, respectively. Even 
without fitting Schechter functions to the data, it is ob- 
vious that the evolution in characteristic luminosity and 
number density between redshifts {z) = 0.6 and {z) = 3.2 



is very moderate if compared to the evolution in the blue 
bands. 



3. Parameterizing the evolution of the LFs 

To better quantify the redshift evolution of the LFs, we 
use the method introduced in FDFLF I. We parameterize 
the evolution of M* and </)* with redshift assuming the 
following simple relations: 



M* (z) = M* + aln(l + z), 
(f)* (z) — (f)'^ {I + z)^ , and 
a (z) — aa = constant. 



(1) 



We then derive the best fitting values for the free param- 
eters a, &, Mq, and (jj^ by minimizing the of 



(2) 



3 = 1 i=l 



4 



for the galaxy number densities in all magnitude and red- 
shift bins simultaneously (for more details see FDFLF I). 
The free parameters of the evolutionary model are con- 
strained for three different cases: 

— Case 1: FDF LFs between redshift (z) ^ 0.65 and 
(z) ^ 2.26 are used. 

— Case 2: FDF LFs between redshift (z) ~ 0.65 and 
(z) ~ 3.21 are used. 

— Case 3: FDF LFs between redshift (z) - 0.65 and 
(z) ^ 2.26 as well as the local LF of iBlanton et alJ 
l|2003() are used. 

As can be seen from F i g. H] t he slopes of the red LFs 
derived by iBlanton et al.l l)2003l SDSS Data release 1) 
are much shallower than those deriv ed in this work and 
a previous w ork of the same author l)Blanton et al-llioOli 
SDSS EDR). IBlanton et all ll2003ll argued that the di ffer- 
ence in the r-band LF (between [Blanton et al.ll200ll and 
IBlanton et"al]l2003|) stems only from the inclusion of lu- 
minosity evohxtion within the covered redshift range. Very 
recentlv iDriver eTal. (2005) showed that the B-band LF 
derived from the Millennium Galaxy Catalogue (MGC, 
which is fully contained within the region of the SDSS) is 
inconsistent with the SDSS z=0 result of Blanton et alJ 
(2003) by more tha n 3g. On the other h and the M* value 
of the B-Band LF of 'Driver et al. (2005) is consistent with 
those derived by iBlanton et al. ( 2001) after that (jf has 
been renormali zed according t o iLiske et al.l l)2003|) . Once 
corrected, the IBlanton' et all (|2aai|) B-band LF agrees 
well with the M GC, the 2dFGRS a nd the ESO Shoe 
Project estimates. iDriver et al.l l|2005() con clude, that the 
discrepancy between the MGC and the IBlanton et alJ 
(|2003) LFs is complex, but predominantly due to a color 
bias within the SDSS. They also conclude, that the color 
selection bias might be a general trend across all filters. 
We compare in Fig. ED the local Schechter functions 





Fig. 6. LFs in the r'-band from low redshift {{z) = 0.3, upper left panel) to high redshift ((z) = 3.2, lower right panel). 
The filled (open) symbols show the LF corrected (uncorrected) for V/Vmax- The fitted Schechter functions for a fixed 
slope a are shown as solid lines. Note that we only fit the LFs from (z) = 0.6 to (z) = 3.2. The parameters of the 
Schechter functions can be found in Table |B| The dotted line represents the local r'-band LF derived from the SDSS 
(Blanton et al. 2003). The Schechter fit for redshift (z) — 0.6 is indicated as dashed line in all panels. 

Table 6. Schechter function fit in the r'-band 



redshift 


interval 


M* (mag) 


4>* (Mpc~ 


') 


a (fixed) 


0.45 - 


- 0.85 


-22.41 +0.23 


-0.18 


0.0025 +0.0002 


-0.0002 


-1.33 


0.85 - 


- 1.31 


-22.67 +0.14 


-0.13 


0.0019 +0.0001 


-0.0001 


-1.33 


1.31 - 


- 1.91 


-22.38 +0.16 


-0.16 


0.0020 +0.0002 


-0.0002 


-1.33 


1.91 - 


- 2.61 


-22.86 +0.13 


-0.11 


0.0019 +0.0002 


-0.0002 


-1.33 


2.61 - 


- 3.81 


-23.06 +0.15 


-0.15 


0.0013 +0.0002 


-0.0001 


-1.33 



as giv en by iBlanton et alJ l(2003f) and iBlanton et al. 

1120011 (f>* has been renormalized according to Liske et al. 
|2003|) for the r', i', and z' bands. Although there is 
a reasonable good agreement between the LFs if one 
focuses on the bright part {M < —20), they disagree at 
fainter magnitudes. On the other hand the slope of the 
LF is strongly dependent on the depth of the survey. 
The flux limit in the r-band selected SDSS survey is 
about rrir < 17.79. A very rough estimate of the absolute 
limiting magnitude at the mean redshift of the survey 
((z) = 0.1) is therefore Mr' « —20. This means that the 
faint-end of the LF as shown in Fig. EH depends on the 
applied completeness correction (see also the discussion 



in lAndreonll20o4> . Therefore, we decide to use only the 
bright part (M < -20) of the SDSS LFs to constrain the 
free evolutionary parameter of Case 3. 

As th e Schechter paramet ers are coupled, and M* 
and (f>* of IBlanton et al ](|m03l) are derived for a different 
slope a, we decide not to use M* and 0* itself, but 
to reconstruct a magnitude binned luminosity function 
out of the S chechter values M*, (/>*, and a given in 
IBlanton et alJ (,,2003.) . Following the method described in 
Sect.0]to visualize the errors of the literature LFs (shaded 
regions in the plots) we derive the 1-magnitude-binned 
LF as shown in Fig. 1101 (solid points). 
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Fig. 7. LFs in the i'-band from low redshift {{z) = 0.3, upper left panel) to high redshift ((z) = 3.2, lower right panel). 
The filled (open) symbols show the LF corrected (uncorrected) for V/Vmax- The fitted Schechter functions for a fixed 
slope a are shown as solid lines. Note that we only fit the LFs from (z) ~ 0.6 to (z) = 3.2. The parameters of the 
Schechter functions can be found in Table |7| The dotted line represents the local i'-band LF derived from the SDSS 
(Blanton ct al. 2003). The Schechter fit for redshift (z) = 0.6 is indicated as dashed line in all panels. 



Table 7. Schechter function fit in the i'-band 



redshift 


interval 


M* (mag) 


(j}* (Mpc~ 


') 


a (fixed) 


0.45 - 


- 0.85 


-22.81 4-0.23 


-0.24 


0.0021 4-0.0002 


-0.0002 


-1.33 


0.85 - 


- 1.31 


-22.91 +0.16 


-0.15 


0.0018 +0.0001 


-0.0001 


-1.33 


1.31 - 


- 1.91 


-22.33 +0.21 


-0.18 


0.0023 +0.0003 


-0.0003 


-1.33 


1.91 - 


- 2.61 


-22.93 +0.14 


-0.13 


0.0019 +0.0002 


-0.0002 


-1.33 


2.61 - 


- 3.81 


-23.06 -1-0.10 


-0.09 


0.0011 +0.0001 


-0.0001 


-1.33 



Fig. EH shows, that a Schechter function fit to the SDSS 
data with a slope of a = —1.33 (as derived from the FDF 
data) results in a reduced xled ^0. This large Xred 
increases the errorbars of the evolutionary parameter 
since we normalize the result of Eq. (O to a Xred ^ ^ 
before calculating the errors. 

The 1(7 and 2a confidence levels of the evolution pa- 
rameters a and b for the different filters and cases are 
shown in Fig. 1111 These contours were derived by project- 
ing the four-dimensional distribution to the a-b plane, 
i.e. for given a and b we use the value of Mq and 0q which 
minimizes the x^(o,^)- For Case 1 (left panel) the er- 



rorbars of a and b are rather large and although the best 
fitting values suggest a redshift evolution we are also com- 
patible (within 2a) with no evolution of M* and 0* . The 
error ellipses for Case 2 (middle panel) are smaller than in 
Case 1 and for the r'-band LF we see a luminosity and a 
density evolution on a 2a level. For the i'-band and z'-band 
LFs we see only a de nsity evolution on a 2 a level. Including 
also the local LF of lBlanton et al.l l)2003(l in the evolution 
analysis as in Case 3 (left panel) we are able to derive 
a and b with higher precision since A/q and 0q are more 
restricted. The luminosity and density evolution is clearly 
visible on more than 2a level. Please note that combining 
different datasets like the FDF and the SDSS can intro- 



Gabasch et al.: The evolution of the luminosity functions in the PDF: II. The red bands 11 




Fig. 8. LFs in the z'-band from low redshift {{z) = 0.3, upper left panel) to high redshift ((z) = 3.2, lower right panel). 
The filled (open) symbols show the LF corrected (uncorrected) for V/Vmax- The fitted Schechter functions for a fixed 
slope a are shown as solid lines. Note that we only fit the LFs from (z) = 0.6 to (z) = 3.2. The parameters of the 
Schechter functions can be found in Table [S] The dotted line represents the local z'-band LF derived from the SDSS 
(Blanton et al. 2003). The Schechter fit for redshift (z) = 0.6 is indicated as dashed line in all panels. 

Table 8. Schechter function fit in the z'-band 



redshift 


interval 


M* (mag) 


(f)* (Mpc~ 


') 


a (fixed) 


0.45 - 


- 0.85 


-23.06 4^0.25 


-0.21 


0.0022 +0.0002 


-0.0002 


-1.33 


0.85 - 


- 1.31 


-23.30 +0.20 


-0.21 


0.0017 +0.0002 


-0.0001 


-1.33 


1.31 - 


- 1.91 


-22.71 +0.18 


-0.17 


0.0020 +0.0003 


-0.0002 


-1.33 


1.91 - 


- 2.61 


-23.19 +0.13 


-0.13 


0.0018 +0.0002 


-0.0002 


-1.33 


2.61 - 


- 3.81 


-23.42 +0.10 


-0.13 


0.0010 +0.0001 


-0.0001 


-1.33 



duce systematic errors due to different selection techniques 
and calibration differences not fully taken into account 
(see also discussion above). Nevertheless, a comparison of 
the FDF LFs with the SDSS Schechter functions in Fig. HI 
shows a relatively good agreement (at the bright end). 
Furthermore, a detailed comparison of the UV LFs of the 
FDF w ith the LF derived in l arge surveys e . g.lWolf et alJ 
l)2003l based on C 0MB0-17), 'Steidel et al .1 lll999l based 
on LBG analvsis). ll^ata et al.. (, 20031 : .Ouchi et all l|2004L 
based o n Subaru Deep F ield /Survey) or pencil beam sur- 
veys e.g'Pol i et al ] lj200l[ based on both HDFs) presented 
in FDFLF I shows good agreement in the overlapping 
magnitude range at all redshifts. We are thus confident. 



that remaining systematic differences (e.g. due to the in- 
fluence of large scale structure; LSS) must be small. 

The values for the free parameters a, b, Mg , and (pQ as 
well as the associated errors can be found in Table El The 
evolution parameters a, 5, Mg , and 4>q derived in Case 1, 
Case 2, and Case 3 agree all within 2(T. Most of the values 
differ only by la or less. 

In Fig. El we illustrate the relative redshift evolu- 
tion of M* for the different filters and different cases, 
whereas the relative redshift evolution of 0* is shown 
in Fig. El Note that a, &, Afg, and 0^ were derived 
by minimizing Eq. (jSJ and not the differences between 
the (best fitting) fines and the data points in Fig. El 
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Fig. 9. Comparison of the red FDF LFs in the redshift range (z) = 0.3 (u pper panels) and {z) = 0.65 (lower pa nels) 
with t he local Schechter functions as derived in the SDSS bv iBlanton et a l. (2003, dotted line), and iBlanton et alJ 
l)2005l dashed line). The filled (open) symbols show the FDF LF corrected (uncorrected) for V/Vmax- 



and Fig. USl As for the blue bands (FDFLF I) the 
simple paramctrization of Eq. ^ is able to describe the 
evolution of the galaxy LFs also in the red bands very well. 

Recently IBlanton et al.l l|2005fl used the data of 
the SDSS Data Release 2 to analyze the very local 
(0.00 < z < 0.05) LF (corrected for surface-brightness in- 
completeness) down to extremely low luminosity galaxies. 
They found, that a Schechter function is an insufficient 
parametrization of the LF as there is an upturn in the 
slope of the LF for Mr — 51og(/iioo) > -18. We there- 
fore compare in Fig. El the red FDF LFs in two redshift 
ranges ((z) = 0.3 and (2) ~ 0.65) w ith the local Schechte r 
func tions as deriv ed in the SDSS bv lBlanton et al.l ((20031) . 
and IBlanton et al. (2005). Considering the small volume 
covered by the FDF in the redshift bin (z) = 0.3 and 
the fact, that we see clustered spectroscopic redshifts at 
z — 0.22, z — 0.33, and z = 0.39, the agreement be- 
tween the LFs and the Schechter functions is relatively 
good for M < —19. For the fainter part, th e measured 
number density disa gree with IBlanton et al.l (jlopS} and 
IBlanton et al.l l|2005|) in all three analyzed bands. If we do 
the same comparison at (z) — 0.65 where the FDF cov- 
ers a relatively large volume minimizing the influence of 



LSS, the measured LFs follow the very local Schechter 
function of Blan ton at a l. (2005.) also in the faint magni- 
tude regime. M oreover, the upturn of the faint-end of the 
LF as fou nd by 'Blant on et all l)2005l) in the SDSS or by 
IPopesso et al. (200 5) in the RASS-SDSS Gala xy Cluster 
Survey (see also IPerez-Gonzalez et al.ll2005l) . is visible 
also in the FDF data (at least at (z) = 0.65). 
This upturn seems to be less pronounced in the UV 
(FDFLF I). A possible reason for this could again be the 
different contribution of the SED-type LFs presented in 
Fig. ^ In the red bands, the difference between the char- 
acteristic luminosities between the LFs for types 1, 2, 3 
and type 4 together with the dominance of the type-4 LF 
at the faint end results in a dip at M ^ —20. 

Although a Schechter functio n is an insu f ficient 
parametrization of the LF derived bv lBlanton et al.l(l200,5l) 
we used their results as local reference point to calculate 
the evolution of the LF in the various bands by minimizing 
Eq. ^ . Due to the upturn of the faint-end of the local LF 
and the fact that our evolutionary model assumes a nor- 
mal Schechter function, the reduced of Eq. ((5)) is of 
the order of 9. As we do not want to increase the number 
of free parameters by using a double Schechter function 
(at higher redshifts the data are not able to constrain a 
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22.30 a = -l,33 
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Fig. 10. The local LFs as given bv lBlanton et alJ lj2n0.'^ for the r' (left panel), i' ( niiddle panel), a nd z' (left panel) 
bands. The dotted lines in all plots represent the best fitting Schechter function of iBlanton et all l|2nn.'^ . The solid 
points and the associated errorbars are derive d by the Schechter v alues and corresponding errors of th e latter (see 
text). The dashed lines represent the result of iBlanton et alJ l)200l|) after renormalizing cj)* according to iLiske et alJ 
l)2003|) . We also fit a Schechter function (solid line) with a fixed slope of -1.33 as derived from the FDF data do the 
LF (solid dots). The corresponding M*, 0*, as well as the reduced of the fit is also given in the figures. 
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Fig. 11. 1(7 and 2ct confidence levels of the parameters a and h in different bands (r', i', and z') for the evolutionary 
model described in the text. Left panel: FDF LFs between redshift (z) ^ 0.65 and {£) ~ 2.26 are used (Case 1). The 
values for a and h can be found in Table Middle panel: FDF LFs between redshift (z) ^ 0.65 and (z) ^ 3.21 are 
used {Case 2). The values for a and b can be fou nd in TablelHl Right panel: FDF LFs between redshift (z) ^ 0.65 and 
(z) ~ 2.26 as well as the local LF of IBlanton et al... (2QIS) are used {Case 3). The values for a and b can be found in 
Table O 



possible upturn in the LF), we increase the errors of a, 6, 
Q. We do this by an appropriate scaling of the 



M*, and 



errors aij of Eq. (|2J) to obtain a reduced of unity. A 
comparison of the evolution parameter a and b with those 
derived in Case 3 shows, that the evolution in the charac- 
teristic luminosity agrees with Case 3, but the evolution 
of the characteristic density decreases from 6 ~ —0.7 to 
b ^ —0.5 being closer to Case 1 and Case 2. However, a 
no-evolution hypothesis can be excluded on a 2g lev el in 
all three bands if the results of IBlanton et al.l ((200,'j ^ are 
used as local reference points. 

If we compare the evolutionary parameters a and 6 
of the red bands with those derived in the blue bands 
(FDFLF I), the following trend can be seen: with increas- 



ing waveband the redshift evolution of M* and (j>* de- 
creases. Furthermore, if we incl ude in our analysis also 
the results obtained in the SDSS l|Blanton et alJl2003(l the 
brightening of M* and the decrease in 0* for increasing 
redshift is still visible in the red bands at more than Su. 



4. Comparison with observational results from 
literature 

To put the FDF results on the evolution of the LFs into 
perspective, we compare them to other surveys using the 
following approach: 

First we convert results from the literature to our cosmol- 
ogy (f^M^O.3, f^A^O.T, and i7o = 70 kms"! Mpc^^). 
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Fig. 12. Relative evolution of M* with redshift. The solid line represent the best fit of the evolutionary model 
according to Eq. Left panel: FDF LFs between redshift (z) ~ 0.65 and (z) ^ 2.26 are used to constrain the 
evolutionary model {Case 1). Middle panel: FDF LFs between redshift (z) ^ 0.65 and (z) ~ 3.21 are used to constrain 
the evoluti onary model (Case 2 ). Right panel: FDF LFs between redshift (z) ~ 0.65 and (z) ^ 2.26 as well as the 
local LF of lBlanton et all ll2003l) ' are used to constrain the evolutionary model {Case 3) (see also text). 




Fig. 13. Relative evolution of 0* with redshift. The solid line represent the best fit of the evolutionary model 
according to Eq. (Q. Left panel: FDF LFs between redshift (z) ~ 0.65 and (z) ^ 2.26 are used to constrain the 
evolutionary model {Case 1). Middle panel: FDF LFs between redshift (z) ^ 0.65 and (z) ^ 3.21 are used to constrain 
the evoluti onary model {Case 2 ). Right panel: FDF LFs between redshift (z) ~ 0.65 and (z) ^ 2.26 as well as the 
local LF of lBlanton et al.l 1 2003^ are used to constrain the evolutionary model {Case 3) (see also text). 



Although this conversion may not be perfect (we can only 
transform number densities and magnitudes but lack the 
knowledge of the individual magnitudes and redshifts of 
the galaxies), the errors introduced in this way are not 
large and the method is suitable for our purpose. Second^ 
in order to avoid uncertainties due to conversion between 
different filter bands, we always convert our data to the 
same band as the survey we want to compare with. Third, 
we try to use the same redshift binning as in the literature. 

To visualize the errors of the literature LFs we per- 
form Monte-Carlo simulations using the AM*, Ai/)*, and 
Aa given in the papers. In cases where not all of these 
values could be found in the paper, this is mentioned 
in the figure caption. We do not take into account any 
correlation between the Schechter parameters and assume 
a Gaussian distribution of the errors AM*, A0*, and Aa. 



From 1000 simulated Schechter functions we derive the 
region where 68.8 % of the realizations lie. The resulting 
region, roughly corresponding to Icr errors, is shaded in 
the figures. The LFs derived in the FDF are also shown as 
filled and open circles. The filled circles are completeness 
corrected whereas the open circles are not corrected. The 
redshift binning used to derive the LF in the FDF as well 
as the literature redshift binning is given in the upper 
part of every figure. Moreover, the limiting magnitude of 
the respective survey is indicated by the low-luminosity 
cut-off of the shaded region in all figures. If the limiting 
magnitude was not explicitly given it was estimated from 
the figures in the literature. 

A comparison of our FDF results with LFs based on 
spectroscopic distance determinations telantgn e^_gjJ 
■2003. .2005.: Tin et al.. .1996. .1997: .Brown et alJ I2OOII: 
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Filter Case a b 
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IShaplev et all I2OOII: lllbert et alJ l2004|) as well as with 
LPs based on photometric redshift s ijWolf et alJ l2003t 
IChen et aLlEoM iDahlen et alJl20n5l) follows: 



\Blanton et all }i200± \200^) : 

In Fig. |5| we compare the red FDF LPs in two redshift 
regimes {{z) = 0.3 and {z) = 0.65) with the local Schechte r 
func tions as derived in the SDSS bv lBlanton et al. I ll200.1) . 
and iBlanton et al.1 l)2005l) . As previously discussed, the 
agreement between the LFs and the Schechter functions 
is relatively good for M < —19. For the fa inter part, the 
measu red num ber density disagr ee with iBlanton et al.l 
1)2003(1 and Blanton et all l|2005|) . If we do the same 
comparison at (z) = 0.65 where the FDF covers a rela- 
tively large volume minimizing the influence of LSS, the 
measured LFs follow the very local Schechter function of 
IBlanton et al. I ^l2005^ also in the faint magnitude regime. 
Note that IBlanton et all l)2005|) explicitly corrected for 
surface-brightness incompleteness when deriving the very 
local LFs. 



\Lin et all \l99A) : 

Despite the small volume covered by the FDF at low 
redshift we comp are in F ig. [HI (left panel) our LF with 
the LF derived bv lLin e~al. (1996 ) in the Las Campanas 
Redshift Survey (LCRS). Their sample contains 18678 
sources selected from CCD photometry in a "hybrid" red 
Kron-Cousins R-band with a mean redshift of (z) '^0.1. 
The solid line in Fig. 1141 (left p a nel) r epresents the 
LF in the R-band from iLin et al.l 1 1996(l whereas the 
filled circles show our V/Vmax corrected LF derived at 
0.15 < z < 0.45. There is a rather large disagreement 
between the LF in the FDF and in the LCRS, which is 
mainly due to the different slope {a = —0.7 for the LCRS) 
but also the FDF galaxy number density at the bright 
end seems to be slightly higher than in the LCRS. This 
might be partly attributed to cosmic variance and/or 



to the selection method. The difference at the faint-end 
is a well known LCRS feature related to their selection 
method which biases LCRS towards early type systems. 
Indeed, our LF for SED type 1 galaxies (tr iangles in 
Fig. I14() shows a very good agreement with iLin et all 

mm) . 



\Lin et all \l99i) : 

Based on 389 field galaxies from the Canadian Network 
for Observational Cosmology cluster redshift survey 
(CNOCl) selected in the Gunn-r-band iLin et all l|l997fl 
derived the LF in the restframe Gunn-r-band. In Fig. [TH 
(right panel) we compare our luminosity function with 
the LF derived bv iLin et al. (199?) in the redshift range 
z = 0.2-0.6. There is a very good agreement between the 
FDF data and the CNOCl survey concerning the LF, 
if we compare only the magnitude range in common to 
both sur-yeys (s haded region). Also the slope derived in 
iLin et all l|l997(l (a = -1.25 ± 0.19, Table 2 of the paper) 
is compatible with the slope in the FDF. 



Brown et al. 




2001) 


Brown et al. 


( 


2001) 



64 deg^ of V and R images to 
measure the local V- and R-band LF. They analyzed 
about 1 250 V fc R selecte d galaxies from the Century 
Survey ('Cell er et al.l Il997l) with a mean spectroscopic 
redshift of (z) ~ 0.06. 

A comparison between the LF of lBrown et all l|200lj) and 
the FDF is shown in Fig. E| for the V-band (left panel) 
and the R-band (right panel). Although the agreement 
is quite good for the bright end, the number density of 
the faint-end is substantially higher in the FDF (while 
the slope of the LF derived in the F DF is a = —1.25, the 
slope derived bv lBrown et"ai] l|200ll) is a = -1.09 ± 0.09 
in the V- as well as in the R-band). 
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Fig. 14. Left panel: Compariso n of the R - band luminosity function of the FDF (filled/open circles, (z) ^ 0.3) with 
the Schechter function derived in Lin et "dl \199 ^.) {{z) - 0.1). The shaded region is based on AM*, A(/)*, and Aa. The 
triangles show the r'-band FDF LF at (z) ^ 0.3 for SE P type 1 g a laxies . Right panel: Comparison of the Gunn-r-band 
LF of the FDF with the Schechter function derived in lLin et al.l ill 99 /I) (z = 0.2-0.6). The shaded region is based on 
AM*, A0*, and Aa. 



Shaplev et al.l Ij200n) analyzed 118 photometrically 
selected LBGs with Kj-band measurements covering an 
area of 30 arcmin^. 63 galaxies have additional J-band 
measurements and 81 galax ies are spect r oscopi cally 
confirmed. Using this sample IShaplev' et all l|200ltl de- 
rived the luminosity function in the restframe V-band 
at redshift of (z) ^ 3.0. Fig. 1161 shows a c o mpari son of 
the V-band LF derived by IShaplev et alJ l)200l|) with 
the LF in the FDF at (z) ^ 3.0. The agreement is very 
good if we again concentrate on the shaded region. On 
the other hand, because of the depth of the FDF we can 
trace the LF two magnitudes deeper and therefore give 
better constraints on the slope of the Schechter function. 
Comparing the faint-end of the FDF LF with the extrap- 
olated Schechter function of IShaplev et al. (2001j) clearly 
shows, that the very steep slope of a = —1.85 is not seen 
in the FDF dataset. 



Ilbert et al\ t200A): 

Rbert'et'lLLl l|2004^ investigated the evolution of the 
galaxy LF from the VIMOS-VLT Deep Survey (VVDS) 
in 5 restframe bands (U, B, V, R, I). They used about 
11000 objects with spectroscopic distance information 
in the magnitude range 17.5 < I < 24.0 to constrain 
the LF to redshift z ~ 2. In Fig. [HI we compare the 
V, R, and I band LF of the FDF with the Schechter 
function derived in the VVDS survey for different redshift 
bins: 0.20 < z < 0.40, 0.40 < z < 0.60, 0.60 < z < 0.80, 



0.80 < z < 1.00, 1.00 < z < 1.30, and 1.30 < z < 2.00. 
Because of the limited sample size of the FDF at low 
re dshift we cou l d not use the same local redshift binning 
as lllbert et alJ <)2004(l . We therefore compare in Fig. El 
(first row) the VVDS Schechter function at (z) ^ 0.3 
(light gray) and (z) ~ 0.5 (dark gray) with the FDF LF 
derived at 0.2 < z < 0.6 as well as in Fig. El (second 
row) the Schechter function at (z) ^ 0.7 (light gray) 
and (z) ~ 0.9 (dark gray) with the FDF LF derived at 
0.6 < z < 1.0. There is a very good agreement between 
the FDF data and the VVDS survey at all redshifts under 
investigation if we compare only the magnitude rang e in 
common to both surveys (shaded region). Ilbert et alJ 
l(2004ll derived the faint-end slope from shallower data 
if compared with the FDF which have only a limited 
sensitivity for the latter. Nevertheless, in all three bands 
the differences between the formal a derived in the FDF 
(ay = —1.25 ± 0.03 and a^'fei' = —1.33 ± 0.03 constant 
in redshift) and in the VVDS are compatible within la 
to 2(T up to redshift z ^ 0.8 (only one bin in the I-band 
LF differs by slightly more than 2a). At higher redshift 
we do not s ee th e steep slope (^ —1.5) as derived by 
lllbert et aD l(2004l) . The circumstance that in the FDF we 
are able to follow the LF about 3-4 magnitudes deeper 
may explain the disagreemen t betw een the extrapolated 
faint-end slope of lllbert et all (|200i) and the FDF result. 



\Wolf etalinnE}: 

In Fig. El we compare the r'-band LF of the FDF 
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Fig. 15. Comparison of the V-band (left panel) and R-band (right panel) LF of the FDF with the local {{z) ^ 0.06) 
Schechter function derived \a \Brown et al\ \200i() . The shaded region is based on AM*, and Aa. 
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Fig. 16. Comparison of the V-band luminosity func- 
tion of the F DF with the Schechter function derived in 
\Shavlev et al. (2001) at {z) ~ 3.0. The shaded region is 
based on AM*, Ac/)*, and Aa, where the cut-off at low lu- 
minosity indicates the limiting magnitude of the sample. 



with the R-band selected luminosity funct ion derived in 
the COMBO-17 survey l)Wolf et alJ boosi) for different 
redshift bins: 0.2 - 0.6, 0.6 - 0.8, 0.8 - 1.0, 1.0 - 1.2. 
Because of the limited sample size of the FDF at low 



re dshift we co u ld not use the same local redshift binning 
as IWolf et all l|200,ll) . We compare therefore in Fig. IH^ 
(upper left panel) the COMBO 17 Schechter function at 
(z) ~ 0.3 (light gray) and (z) ~ 0.5 (dark gray) with the 
FDF LF derived at 0.2 < z < 0.6. There is a very good 
agreement between the FDF data and the COMBO-17 
survey at all redshifts under investigations if we compare 
only the magnitude range in common to both surveys. 
Although the number density of the FDF seems to be 
shghtly higher for the restframe UV LF (FDFLF I), this 
is not the case i f we compare the LF in the R-band. 
IWolf et 3,1.1 (|2Q03) derived the faint-end slope from rela- 
tively shallow data which have only a limited sensitivity 
for the latter. This may explain the disagreement bet ween 
the extrapolated faint-end slope of I Wolf et al.l l)2003|) and 
the FDF result. 



\r,hen et al\ fEoM): 

The galaxy sample analyzed by IChen et al.l l)2m);il ) con- 
tains ^ 6700 H-band selected galaxies (within 847 
arcmin^) in the HDFS region with complementary op- 
tical U, B, V, R, and I colors, and ^ 7400 H-band 
selected galaxies (within 561 arcmin^) in the Chandra 
Deep Field South region with complementary optical V, 
R, I, and z' colors. The galaxy sample is part of the 
Las Campanas Infrared Survey (LCIR Marzk e et al1ll999t 
iMcCarthv eTal] 1200 ll) and based on photometric red- 
shifts. 

Fig. 1191 shows a c omparison o f the R -band luminosity 
function derived by IChen et alJ l)2003|) with the LF in 
the FDF for three different redshift bins: 0.50-0.75 (left 
panel), 0.75-1.00 (middle panel), and 1.00-1.50 (right 
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Fig. 17. Comparison of the V (left pane l s), R (middle panels), and I (right panels) band LF of the FDF with the 
Schechter function derived in lllbert et all \200J ) (VVDS) at 0.20 < z < 0.40 (first row, Ught gray), 0.40 < z < 0.60 
(first row, dark gray), 0.60 < z < 0.80 (second row, light gray), 0.80 < z < 1.00 (second row, dark gray), 1.00 < z < 1.30 
(third row), and 1.30 < z < 2.00 (fourth row). The shaded regions of all plots with z < 1 are based on AM*, A0*, 
and Aa. Only in the two high redshift bins (third and fourth row) the shaded regioi i is based on l y on A M* and A(j>* . 
Please note that we use the average error if the upper and lower values reported bv lllbert et alJ |20^ disagree. 
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Fig. 18. Comparison of the LF in the r'-band of the FDF with the Schechter function derived in l TyoZf ~ 
0.2 < z < 0.4 (upper left panel, light gray), 0.4 < z < 0.6 (upper left panel, dark grey), 0.6 < z < 0.8 (upper right 
panel), 0.8 < z < 1.0 (lower left panel), 1.0 < z < 1.2 (lower right panel). The shaded regions of nearly all plots are 
based on AM*, A0*, and Aa. Only in the highest redshift bin (lower right panel) the shaded region is based only on 
AM* and A^*. The cut-off at low luminosity indicates the limiting magnitude of the sample. 



panel). There is a good agreement be tween the FDF LF 
and the Schechter function derived bv lChen et alJ l)2003|) 
in the lowest redshift bin (z = 0.50-0.75) if we compare 
only the magnitude range in common to both surveys. At 
intermediate redshift (z = 0.75-1.00) the number density 
of the bright e nd of the FDF LF is slightly higher than 
in lChen et al.l l|2003|) . On the other hand, for the highest 
redshift bin (z = 1.00 -1.50) the number density of the 
bright end derived by IChen et all l|200,lf) roughly agrees 
with the results obtained in the FDF. 



Dahlen et al i200A} : 

Dahlen et all l)2005l) used HST and ground-based U 
through Kg photometry in the GOODS-S Field to mea- 
sure the evolution of the R-Band luminosity function out 
to z 2. They combine a wider area (^ 1100 arcmin^), 
optically selected {R < 24.5) catalog with a smaller area 
(~ 130 arcmin^) but deep NIR selected {K^ < 23.2) 
catalog. Distances are based on photometric redshifts 
with an accuracy of Az/{zspec -I- 1) ~ 0.12 0.06 after 
excluding ~ 3% of outliers). To determine the restframe 
R-band galaxy luminosity function out to z ~ 2 they 
used the deep Kg selected catalog. A comparison between 




Fig. 19. Compari son of the luminosity function in the R-band of the FDF with the Schechter function derived in 
\Chen et ol\ \200A) : 0.50 <z< 0.75 (left panel), 0.75 <z< 1.00 (middle panel), and 1.00 <z< 1.50 (right panel). 
The shaded region is based on AM*, A4>* and Aa for 0.50 < z < 0.75 (left panel). For 0.75 < z < 1.00 (middle panel) 
and 1.00 < z < 1.50 (right panel) the shaded region is based only on AM* and A0*. The cut-off at low luminosity 
indicates the limiting magnitude of the sample. 



the R-band LF of iDahlen et all l)2005|) and the FDF 
is shown in Fig. 1201 There is a very good agreement 
in nearly all redshift bins. Only at 0.82 < z < 1.12 and 
1.37 < z < 1.59 the characteristic density in the FDF 
seems to be slightly higher. 

To summarize we can say, that the LFs derived in the 
FDF in general show a very good agreement with other ob- 
servational datasets from the literature. At the bright end 
of the LF most of the datasets agree within la. Differences 
between the extrapolated Schechter function of the liter- 
ature and the measured faint-end in the FDF can be at- 
tributed to the shallower limiting magnitudes of the other 
surveys. 

5. Comparison with model predictions 

As discussed for example in lBenson et al.l l)2003j) different 
physical processes are involved in shaping the bright and 
the faint-end of the galaxy LF. Therefore, it is interesting 
to compare LFs predicted by models with observational 
results to better constrain those processes. In this section 
we compare the R-band and I-band LFs in different 
redshi ft bins with model predictions of iKauffmann et'al] 

In Fig. 12] we show the R-band luminosity function of the 
FD F together with the sem i-analytic model predictions 
by IKauffmann et all \im ^^ for {z) ~ 0.20, {z) ~ 0.62, 
(z) ~ 1.05, whereas in Fig. [53 we show the I-band LF 
in the redshift bins (z) - 0.20, (z) - 0.62, (z) - 1.05, 
(z) ~ 1.46, (z) ~ 2.12, and (z) ~ 2.97. For the R-band no 
semi-analytic model predictions are available for redshifts 
larger than (z) ^ 1.05. 

There is a good agreement between model (solid lines) 

^ The models were taken from: 
jhttp:/ /www.mpa-garching.mpg.de/Virgo/data_download.htnil| 



and measured LFs in the R-band. Also for the I-band 
there is a good agreement between the models and the 
luminosity functions derived in the FDF up to redshift 
(z) ~ 1.46 (of course at z « the model is tuned to repro- 
duce the data). At (z) > 1.46 the discrepancy increases 
as the model does not contain enough bright galaxies. 
Unfortunately, the models only predict luminosities for 
massive galaxies and because of lack of resolution do not 
predict galaxy number densities for faint galaxies. 



6. Summary and conclusions 

In this paper we use a sample of about 5600 I-band se- 
lected galaxies in the FORS Deep Field down to a lim- 
iting magnitude of I = 26.8 mag to analyze the evolu- 
tion of the LFs in the r', i', and z' bands over the red- 
shift range 0.5 < z < 3.5, thus extending the results pre- 
sented in FDFLF I to longer wavelengths. All the results 
are based on the same catalog and the same state of 
the art photometric redshifts (Az/(z5pej, -I- 1) < 0.03 with 
only ~ 1% outliers) as in FDFLF I. The error budget of 
the luminosity functions includes the photometric redshift 
error of each single galaxy as well as the Poissonian error. 
Because of the depth of the FDF we can trace the LFs 
deeper than most other surveys and therefore obtain good 
constraints on the faint-end slope a of the LF. A de- 
tailed analysis of a leads to similar conclusions as found in 
FDFLF I for the blue regime: the faint-end of the red LFs 
does not show a large redshift evolution over the redshift 
range 0.5 < z < 2.0 and is compatible within la to 2a with 
a constant slope in most of the redshift bins. Moreover, the 
slopes in r', i', and z' are very similar with a best fitting 
slope of a = —1.33 ± 0.03 for the combined bands and 
redshift intervals considered here. 

Interestingly, an analysis of the slope of the LFs as a func- 
tion of wavelength shows a prominent trend of a to steepen 
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Fig. 20. Comparison of the R-band LF of the FDF with the Schechter function derived m lPahlen et al\ }200A) at 
0.10 <z< 0.82, 0.82 <z< 1.12, 1.12 <z< 1.37, 1.37 <z< 1.59, 1.59 <z< 1.80, and 1.80 <z< 2.00 (from upper 
left to lower right panel). The shaded region is based on AM*, A0* and Aa for all panels. The cut-off of the shaded 
region at low luminosity indicates the limiting magnitude of the sample. 



with increasing wavelength: a[/y&„' = —1.07 ± 0.04 
ag,&s = -1.25 ±0.03 = -1.33 ± 0.03. To 

better understand this wavelength-dependence of the LF 
slope, we analyze the contribution of different galaxy types 
to the overall LF by subdividing our galaxy sample into 4 
typical SED types with restframe U-V colors between 2.3 
- 1.9, 2.0 - 1.6, 1.6 - 0.9, and 0.9 - for SED type 1, 2, 
3, and 4, respectively. Therefore, in a rough classification 
one can refer to SED types 1 and 2 (SED type 3 and 4) 
as red (blue) galaxies. 

Although in the UV regime the overall LF is completely 
dominated by extremely star-forming galaxies, the overall 
LF in the red regime is mainly dominated by early to late 
type galaxies at the bright end, but extremely star-forming 
galaxies at the faint-end. The relative contribution of the 
different SED type LF to the overall LF clearly changes as 
a function of analyzed waveband resulting (at the depth 
of the FDF) in a steeper slope for the overall LF in the 
red regime if compared to the blue regime. 
To quantify the contribution of the different SED types 
to the total luminosity density, we derive and analyze the 
latter in the UV and in the red bands as a function of 
redshift: The contribution of type 1 and 2 galaxies to the 



UV LD is negligible at all analyzed redshifts as SED type 
3 and 4 completely dominate. On the other hand, the rela- 
tive contribution to the overall luminosity density of type 
1 and 2 galaxies is of the same order or even exceeds those 
of type 3 and 4 in the red bands. 

We investigate the evolution of M* and </>* (for a fixed 
slope a) by means of the redshift parametrization intro- 
duced in FDFLF I. Based on the FDF data [Case 1 and 
Case 2), we find only a mild brightening of M* and de- 
crease of 0* with increasing redshifts in all three analyzed 
wavebands. If we follow the evolution of the characteristic 
luminosity from (z) ~ 0.5 to (z) ^ 3, we find an increase 
of ~ 0.8 magnitudes in the r', and '--^ 0.4 magnitudes in 
the i' and z' bands. Simultaneously the characteristic den- 
sity decreases by about 40 % in all analyzed wavebands. 
We compare the LFs with previous observational datasets 
and discuss discrepancies. As for the blue bands, we find 
good/ very good agreement with most of the datasets espe- 
cially at the bright end. Differences in the faint-end slope 
in most cases can be attributed to the shallower limiting 
magnitudes of the other surveys. 

We also compare our results with predictions of semi- 
analytical models at various redshifts. The semi-analytical 



22 



Gabasch et al.: The evolution of the luminosity functions in the FDF: II. The red bands 



M 10-3 



10-5 



— Kauffmann et al, 19 
o • fDf < z g 0.4 




-24 



-22 



-20 



M 10-3 



10-5 



Kauffmann et al, 1999 <z> -0 62 
FDf 0,-1 < z g 0.8 




-24 -22 



-20 



I 

-IE 



M 10-3 

E 



10-' 



— Kauffmann et al. 1999 <z> -1.05 
o • fDl' 0.8 < z g 1.2 




-24 



-22 



-20 



Fig. 21. Comparison of the R-band LF of the FDF with predictions based OTi \Kauffm,ann et 

!ID \M2^ (solid line): 

(z) ^ 0.20, (z) ^ 0.62, and (z) 1.05, (from left to right panel). The filled (open) symbols show the LF corrected 
(uncorrected) for V/Vmax- The drops of the t heoretical curves towards the faint-end is caused by the limited mass 
resolution of the models, see lKauffmann et al. for details. 



models predict LFs which describe the data at low redshift 
very well, but as for the blue bands, they show growing 
disagreement with increasing rcdshifts. Unfortunately, the 
models only predict luminosities for massive galaxies and 
therefore, a comparison between the predicted and ob- 
served galaxy number densities for low luminosity galaxies 
{L <L*) could not be done. 
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